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We report on neutron diffraction, thermal expansion, magnetostriction, dielectric, and specific heat 
measurements on polycrystalline FeCr 2 S 4 in external magnetic fields. The ferrimagnetic ordering 
temperatures T c ~ 170 K and the transition at T 00 ~ 10 K, which has been associated with orbital 
ordering, are only weakly shifted in magnetic fields up to 9 T. The cubic lattice parameter is found to 
decrease when entering the state below T 00 . The magnetic moments of the Cr- and Fe-ions are reduced from 
the spin-only values throughout the magnetically ordered regime, but approach the spin-only values for 
fields >5.5 T. Thermal expansion in magnetic fields and magnetostriction experiments indicate a 
contraction of the sample below about 60 K. Below T QO this contraction is followed by a moderate 
expansion of the sample for fields larger than ~4.5 T. The transition at T 00 is accompanied by an anomaly 
in the dielectric constant. The dielectric constant depends on both the strength and orientation of the 
external magnetic field with respect to the applied electric field for T < Too A linear correlation of the 
magnetic-field- induced change of the dielectric constant and the magnetic-field dependent magnetization is 
observed. This behaviour is consistent with the existence of a ferroelectric polarization and a multiferroic 
ground state below 10 K. 

The interplay of orbital, lattice, and spin degrees of freedom is of fundamental importance in understanding 
the enormous variety of phenomena and ground states in transition-metal compounds 1 3 . In particular, the 
effects of short and long-range orbital ordering or of dynamic and static cooperative Jahn-Teller (JT) 
distortions on the magnetic ground state have been in the focus of condensed-matter physics throughout the 
research activities on, for example, cuprates and manganites. The respective JT active magnetic ions Cu 2+ and 
Mn 3+ usually exhibit a strong JT coupling in octahedral environment and contribute to cooperative JT distortions 
and orbital ordering far above room temperature, e.g. in KCuF 3 and LaMn0 3 4 " 6 . Similar electronic configurations 
are found in the case of compounds containing Fe 2+ or Cr ,+ -ions in a tetrahedral crystal-field 2,71 ", but the 
electron-lattice coupling is generally weaker and competes with spin-orbit coupling and exchange interactions, 
leading to strong fluctuations and frustration effects. 

FeCr 2 S 4 is a sulfur spinel that has been reported previously to exhibit a complex interplay of orbital correlations 
and magnetism, including ferrimagnetic ordering of the Cr and Fe sublattices below T c ~ 170 K 11,12 and a 
transition at T 00 = 10 K, which has been assigned to originate from orbital ordering 7,12 " 19 . In addition, an 
anomaly in the low-field magnetization at T M = 60 K revealed a change in the magnetic configuration 12,20 " 25 , 
which has recently been attributed to the formation of a non-collinear (possibly helical) spin configuration with 
an incommensurate modulation involving three different Fe sites as indicated by ^iSR, Mossbauer, and time- 
resolved magneto-optical Kerr effect measurements 26 2S . The same temperature has also been associated with the 
onset of short-range orbital order (orbital liquid) and dynamic JT distortions, suggesting a mutual influence of 
spin configuration and orbital correlations 12 . 

At room temperature FeCr 2 S 4 crystallizes in the spinel structure with space group Fd3m with eight formula 
units per unit cell. The Cr 3+ -ions (3d 3 , S = 3/2) occupy the octahedral sites and the Fe 2+ -ions (3d 6 , S = 2) are at the 
centre of S 2 ~ tetrahedron ions, forming a diamond lattice of two interlacing FCC sublattices 2 '. While the Cr 
sublattice is dominated by ferromagnetic exchange which is established through the 90° Cr-S-Cr bond, the 
exchange coupling within the Fe sublattice is rather weak and overruled by a much stronger antiferromagnetic 
coupling to the Cr-ions, resulting in the ferrimagnetic arrangement. As the Cr 3+ -ions with a half-filled t^ground- 
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Figure 1 | Results of the neutron diffraction temperature dependency study, (a) Rietveld refinement of high resolution data collected at 4 K indicates a 
cubic spinel crystal structure and collinear antiparallel spin alignment, (b) High intensity neutron diffraction measured from T = 4 K to 300 K 
represented in a contour plot with a logscale color map (counts per 10 min) to enhance subtle modifications, revealing T c ~ 180 K. (c) Temperature 
dependency of the refined cubic lattice parameter. The dashed line corresponds to the behavior of a normal anharmonic solid, from which FeCr 2 S 4 
deviates due to short range orbital order below Tm ~ 60 K. Inset: Contraction of the lattice below Too = 10 K. (d) Temperature dependency of the 
magnetic moment parameter of the Fe- and Cr-ions. Above 200 K the magnetic phase was not included in the refinement. In (c) and (d) the solid lines 
correspond to data collected on Wombat, and black markers correspond to high resolution data collected on the instrument Echinda. 



state configuration do not possess orbital degrees of freedom, the 
orbital physics in this compound is dominated by the Fe 2+ -ions. 
Because of the internal exchange field below T c , the five-fold spin 
degeneracy of the lower-lying^ states is lifted and the level splits into 
five orbital doublets. The lowest level is coupled to lattice vibrations, 
inducing a dynamic JT effect before an eventual formation of a static 
cooperative JT distortion of the sulfur tetrahedra below T = 10 K 29,30 . 
Previous specific heat measurements show a A -like anomaly at T 00 , 
which was interpreted as evidence for an orbital order at low tem- 
perature 12 . However, an unambiguous signature of a structural JT 
distortion, predicated by a deviation from cubic symmetry, remains 
evasive even in high-resolution structural x-ray studies. Thermal 
expansion measurements have revealed a contraction of the lattice 
at T 00 12 . 

The influence of an external magnetic field on the magnetization 
of FeCr 2 S 4 in the orbitally ordered phase was reported recently and it 
was proposed that there exists a new magnetic phase for fields larger 
than 5.5 T, which might be related to the restoration of a commen- 
surate collinear spin structure 31 . Related compounds FeCr 2 0 4 and 
CoCr 2 0 4 (with orbital degrees of freedom) have been found to be 
multiferroic 32,33 , and CoCr 2 0 4 exhibits unconventional properties as 
a function of applied magnetic field. Therefore, we focus in this study 
on the influence of external magnetic fields on the structural, mag- 
netic, and dielectric properties of FeCr 2 S 4 in the orbitally ordered 
phase. We identify a decrease of the cubic lattice parameter below 
T 00 anQ a step-like increase of the magnetic moments of the Cr-and 
Fe-ions upon application of an external magnetic field. Moreover, 
our dielectric measurements in magnetic fields provide evidence for 
the existence of a ferroelectric polarization in FeCr 2 S 4 . All investi- 
gated quantities show changes in magnetic fields of about 5 T. 



Results and discussion 

Neutron diffraction. High resolution neutron diffraction patterns 
were collected on the instrument Echidna, located at the Bragg 
Institute, ANSTO, at temperatures of T = 4 K, 15 K, 70 K, and 
100 K. Figure 1(a) presents raw data collected at 4 K, together 
with the final Rietveld refinement results. Refinements of the 
diffraction patterns were performed in the high- symmetric cubic 
spinel structure with the space group Fd3m (no. 227). The atoms 
are located at the following Wyckoff positions: Fe at 8a (1/8, 1/8, 
1/8), Cr at 16d (1/2, 1/2, 1/2), and S at 32e (x, x, x). Instrumental 
parameters were fitted using the diffraction pattern measured at 
4 K, and fixed for the refinement of the succeeding datasets. Fitted 
structural parameters of the Rietveld refinements for 4 K and 
15 K are listed in Table I. Within the resolution limits of the 
instrument the sulfur position parameter x did not change 
between 4 and 15 K, but there is a notable decrease of the cubic 
lattice parameter between T = 4 K and 15 K as shown in the inset 
of Fig. 1(c). 

Therefore, we further investigated this trend in a detailed neutron 
diffraction study. Diffraction patterns were collected on the high 
intensity instrument WOMBAT, located at the Bragg Institute, 
ANSTO, at temperatures from 4 K to 300 K in 1 K steps, which 
revealed the transition to ferrimagnetic order below T c ~ 180 K. 
The full 4 K to 300 K diffraction dataset was systematically refined to 
the cubic spinel space group Fd3m. Initial crystal and magnetic phase 
parameters were set using the results of the refinement of the high- 
resolution data. Instrumental characteristics such as peak shape 
parameters were fitted for the 4 K dataset and fixed for succeeding 
temperature refinements. Above 200 K the magnetic phase was not 
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included in the refinement. The refined cubic lattice parameter and 
magnetic moments of the Cr- and Fe-ions are plotted in Fig. 1(c) and 
1(d), respectively, which demonstrate the consistency between the 
temperature dependency of the high-intensity and high-resolution 
neutron diffraction studies. The insert of Fig. 1(c) reveals a distinct 
reduction of the cubic lattice parameter below 10 K. This contraction 
can be interpreted as a microscopic signature of the antiferro-orbital 
long-range order below T QO = 10 K. In Fig. 1(c) we show that the 
temperature-dependent trend in the lattice parameter below T c to 
T M at 60 K is well described by normal anharmonic lattice dynamics, 
and was fit accordingly to V(T) = V Q [1 + A/[exp(0/T) — 1]], where 
V 0 represents the cell volume extrapolated to 0 K, 0 is an averaged 
Debye temperature (0 = 481 K in the ferrimagnetic region), and A is 
a fitting constant. In the transition from paramagnetic to ferrimag- 
netic order at T c a deviation in the lattice parameter contraction is 
found to form as a result of magnetostriction. More intriguing is a 
deviation from the fitted curve for a normal anharmonic solid below 
60 K, which indicates enhanced shrinking of the lattice parameter. 
To our knowledge these deviations in structural characteristics have 
not been microscopically identified using neutron diffraction prior to 
this study. However, they are consistent with previous bulk tech- 
niques such as thermal expansion measurements 12 . The behavior 
below T M ~ 60 K has previously been associated with short-range 
orbital fluctuations 12 and the emergence of a non-collinear spin 
configuration 26 . 

The temperature dependency of the magnetic moment of the Fe- 
and Cr-ions as determined from the refinement of the neutron dif- 
fraction data are presented in Fig. 1(d). The magnetic phase 
symmetry was defined to align the Fe and Cr moments antiparallel 
along the structural c-direction as our data indicates that the spins are 
aligned along the cubic crystal axes. No additional Bragg reflections, 
which would indicate a weak spin canting or incommensurate struc- 
ture, were detected within the experimental resolution. Simulations 
of a possible spin canting did indicate that additional transversal 
magnetic components of no larger than 0.3 /(s/Fe 2+ -ion could exist 
below 60 K. The obtained values of — 3.69 fi B and 2.64 ji B for the Fe- 
and Cr-ions at 4 K, respectively, are consistent with values recently 
determined in a neutron diffraction study 18 , along with other experi- 
mental and theoretical approaches, which are compared in Table II. 
In particular, density-functional based calculations found reduced 



spin moments of —3.3 ji B and 2.7 fi B for both types of ions and an 
orbital moment of —0.13 [i B on the Fe-ions 37,38 . This is in agreement 
with our resulting total magnetic moments within experimental 
error. The orbital contribution of the Fe-ions is argued to originate 
from spin-orbit coupling 38 , which is in agreement with the single-ion 
picture, where the second-order contribution of spin-orbit coupling 
yields an effective g- factor ofg = 2.1 39 . The successful refinement of 
the FeCr 2 S 4 magnetic structure indicated that the spin-only 
moments are reduced from expected values of —4.2 n% for Fe 2+ with 
S = 2 andg = 2.1 and 3.0 fi B for Cr 3+ with spin S = 3/2 andg = 2. 
Moreover, recent reports have indicated modifications to the bulk 
magnetic properties of FeCr 2 S 4 under large applied magnetic fields in 
this temperature range 31,40 . 

Neutron diffraction was performed on FeCr 2 S 4 under external 
applied magnetic fields of up to 7.5 T in order to investigate a poten- 
tial link between orbital order, the magnetic and crystal structures 
and external fields. In particular, deviations in the magnetic and 
crystal structure were studied by performing Rietveld refinement 
of the lattice and magnetic moment parameters of the Fe 2+ and 
Cr 3+ -ions. Neutron diffraction patterns were collected on the instru- 
ment Wombat at ANSTO under applied magnetic fields from 3 T up 
to 7.5 T in 0.5 T steps at temperatures from 5Kuptol5KinlK 
steps. For the Rietveld refinements, the collected diffraction patterns 
were treated using the cubic spinel structure and collinear antipar- 
allel magnetic structure as determined previously, since there was no 
evidence for additional Bragg peaks within the experimental resolu- 
tion. Crystal and magnetic structural parameters were initially set 
based upon the results of the zero-field refinement results (see 
Table I). Under an initial applied field of 0.5 T, the Fe- and Cr-ion 
spins structure of cubic FeCr 2 S 4 preferentially rotates to align along 
the crystal axis in the applied magnetic field direction. Additional 
parameters to account for this preferential orientation of the mag- 
netic phase were implemented in the refinement, resulting in a sig- 
nificant improvement to the fitting residuals. 

The obtained lattice parameter and magnetic moments are pre- 
sented in Fig. 2. In Fig. 2(a) the reduction of the cubic lattice para- 
meter below Too observed in zero-field remains visible under 
applied magnetic fields up to 7.5 T. This indicates that the phase 
below Too remains stable with respect to even relatively large 
external magnetic fields. A clear increase in the obtained magnetic 
moment of the Fe 2+ and Cr 3+ -ions occurs at around 4.5 T 
(Fig. 2b,c,d), with spin-only values of —4.2 fi B and 3.0 f.i B 
approached for fields above 5.5 T, respectively. 

The reduced magnetic moments observed in neutron diffraction 
in zero-field can be explained due to the existence of a transverse 
quasi-paramagnetic component of the magnetic moment of the 
Fe 2+ -ions. According to previous publications, they are a con- 
sequence of the orbital moment of the Fe-ions which causes an 
magnetocrystalline anisotropy with an estimated magnitude of about 
2-4 T in the magnetically ordered regime T < T C 2SJ7JS . A possible 
scenario to understand the observed increase is that these transverse 
components are suppressed in the external magnetic field and the 
magnetic moments reach the full spin-only values of the longitudinal 
magnetization. The strength of the necessary magnetic field is in 
good agreement with reports of a change in slope of the magnetiza- 
tion as a function of field 31,40 and with the order of magnitude of the 



Table II | Comparison of the derived magnetic moments of the Fe- and Cr-ions and the total magnetic moment per formula unit (in units of the 

Bohr magneton /<b) from neutron diffraction (ND), x-ray spectroscopy (XPS), magnetization measurements (Mag.), and theoretical 
approaches (Calc.) for lowest temperatures 

This Work ND[18] ND[11] XPS[34] Mag.[35] Mag.[21] Mag.[31] Calc.[36] Calc.[37] 

Fe 2+ -3.69 -3.52 -4.2 -3.89 - - - -3.70 -3.40 

Cr 3+ 2.64 2.72 2.9 2.71 - - - 3.11 2.66 

Total/f.u. 1.59 1.92 1.6 1.53 1.52 1.84 1.6 2.0 1.92 
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Figure 2 | Results of the magnetic field dependency study have been plotted as contour maps, (a) The cubic lattice parameter shows a clear contraction at 
T = 10 K for all magnetic fields studied. The refined (b) bulk magnetic moment, (c) Fe-ion moment and (d) Cr-ion moment parameters reveal 
a clear increase at around 4.5 T, which is interpreted as the realization of full spin-only values under an applied magnetic field of the order of the 
anisotropy field, as described in the text. 



field necessary to overcome the magnetocrystalline anisotropy in 
FeCr 2 S 4 28 . Therefore, we believe that for fields larger than 5.5 T 
FeCr 2 S 4 becomes an isotropic ferrimagnet. At the same magnetic 
field strength the magnetostriction and the dielectric properties also 
show distinct changes, as discussed in the following. 

Thermal expansion, magnetostriction, and specific heat. The 

temperature dependence of the thermal expansion (L(T) — L 2 ook)/ 
i2ooK °f me sample's thickness L in the direction of the applied 
magnetic fields of a dense polycrystalline sample obtained by 
spark-plasma sintering are shown in Fig. 3(a) together with zero- 
field data reported in Ref. 12. In the paramagnetic phase at 200 K we 
do not expect any significant magneto-strictive effects enabling us to 
compare directly the effect of the external magnetic field, by 
normalizing the data to the thickness L 2 ook of the sample at 200 K. 
The temperature dependence in zero field clearly shows a 
contraction of the sample upon entering the state below 10 K 12 , 
which is in agreement with the temperature evolution of the lattice 
parameter in Fig. 1(c). 

The effect of the external magnetic field above 100 K is small, but a 
clear broadening of the anomaly at the ferrimagnetic transition is 
found in the linear thermal-expansion coefficient a L = \IL{dLldT) 
for the direction parallel to the applied magnetic field (Fig. 3(b)). In 
the temperature range 60 £ T £ 130 K only slight deviations from 
the zero-field curve can be detected. However, below about 60 K the 
sample contracts much stronger when a magnetic field is applied. 
This temperature coincides with the temperature T M where a trans- 
ition from a collinear to a non-collinear spin configuration with three 
different Fe sites has been reported by iiSR and where orbital fluctua- 
tions are assumed to set in 12,26,27 . The relative change of the thermal 
expansion upon entering the ground state below 10 K is similar to 



the zero-field data, but under applied fields of 6 T and 12 T the 
sample expands again towards lowest temperatures. This is clearly 
seen in the negative thermal expansion coefficient for the transition 
in magnetic fields. The temperatures of the minima in the thermal 
expansion in magnetic fields are used as a measure of the orbital- 
ordering transition temperatures. 

The relative magnetostriction effect of the sample (1(H) — L H=0 )/ 
Lh=o) is shown in Fig. 3(c) for selected temperatures, which are also 
indicated by arrows in the inset of Fig. 3(a). The values at 12 T are 
compared with the values (large solid symbols in Fig. 3(c)) obtained 
by subtracting the corresponding temperature dependencies in 
Fig. 3(a). The different data sets agree well. The deviations at lower 
temperatures might indicate small variations in temperature during 
the different measurements. At high temperatures the field effects are 
small and the almost field independent zero-magnetostriction curve 
at 70 K indicates a sign change as expected from the temperature 
dependence of the thermal expansion. Accordingly, the largest con- 
traction is observed at 12 K, just above the orbital ordering trans- 
ition. At both 12 K and 25 K the lattice contracts strongest in field 
from zero to about 2 T, undergoes a very broad shallow minimum, 
and tends to saturate up to the highest measured fields. In the orbi- 
tally ordered state a broad but clear minimum develops at around 3- 
4 T, signaling the change from a contraction to an expansion of the 
sample in higher magnetic fields. The external fields B cl and B c2 
associated with this change are estimated by the values where the 
magnetostriction curve increases visibly from the minimal value 
(indicated as a constant solid line) and by the inflection points, 
respectively. 

From the neutron scattering results in magnetic field, it is clear 
that the volume contraction of the lattice below T OQ is still present in 
magnetic fields up to 7.5 T. The strong effects seen in the mac- 
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Figure 3 | Temperature dependency of (a) the thermal expansion normalized to the value at 200 K and (b) the thermal expansion coefficient a in 
magnetic fields of 0, 6, and 12 T. Inset in (a): thermal expansion for T £ 80 K on an enlarged scale. Arrows indicate selected temperatures where the 
differences in the thermal expansion in constant magnetic fields are compared to the thermal expansion vs. magnetic field normalized to the zero-field 
value measured at different temperatures shown in (c). The values derived from the temperature dependencies as described in the text are indicated 
by large spheres. The arrows at H cl and H c2 mark the fields associated with the beginning expansion along the direction of the magnetic field, 
(d) Temperature dependency of the specific-heat in zero and 9 T plotted as C P /Tvs. T showing the ferrimagnetic and orbital ordering transitions 
at T c and T 00 , respectively. Inset: Changes of the anomaly at the orbital-ordering transition for 0, 1, 6, and 9 T. Solid lines indicate linear fits to estimate 
the transition temperatures as described in the text. 



roscopic thermal expansion measurement must therefore be related 
to the existence of magnetic domains in the investigated sample. In 
fact, the largest changes occur in fields smaller than 2 T (see Fig. 3)(c) 
at temperatures below T M = 60 K, where strong domain formation 
has been reported 22 in agreement with the observation of small hys- 
teretic effects. Within the magnetic field range of our study, the 
contraction of the sample seems to saturate above 3-4 T and above 
Tqo> indicating that the magnetic domain structure is stabilized. In 
the orbitally ordered state, however, an expansion of the sample sets 
in between B cl and B c2 . As mentioned above this field range coincides 
with the increase of the magnetic moments of the Fe- and Cr-ions 
and reflects the order of magnitude of the magnetocrystalline 
anisotropy. 

In Fig. 3(d) we compare the temperature dependence of the spe- 
cific heat in the representation C P /T vs. T in zero and 9 T fields. 
While the sharp anomaly at the ferrimagnetic transition at 165 K 
in zero field is strongly broadened at 9 T, the changes of the 1-like 
anomaly around 10 K at the orbital- ordering transition (shown in 
the inset for several external magnetic fields) exhibits only a small 
shift and a transfer of associated entropy to higher temperatures with 
increasing magnetic field. This indicates that the orbital configura- 
tion of the ground state is rather stable with respect to the applied 
magnetic field, but critical fluctuations are favored to occur at higher 
temperatures already. A shift of the anomaly and an entropy transfer 
at 14 T has been reported 40 , which is similar to the one we observe at 
9 T. The temperatures of the orbital-ordering transition at 9.2, 9.5, 



9.9, and 10.3 K for 0, 1,6, and 9 T, respectively, were evaluated from 
the Cp/T dependencies in terms of the temperatures where the data 
starts to deviate from a linear fit at temperatures above the maximum 
(see inset of Fig. 3(d)). The ferrimagnetic transition temperature in 
zero-field is determined by the same procedure as T c = 165 K. At 
9 T the transition temperature has been evaluated as the intersection 
point of two linear extrapolations below and above the broad max- 
imum (see Fig. 3(d)). 

Dielectric spectroscopy. In this section we present detailed measure- 
ments of the dielectric constant as a function of temperature and 
external magnetic fields specifically focusing on the occurrence of 
magnetoelectric effects. Temperature and frequency dependent of 
the dielectric loss has been reported earlier. In single crystals where 
the transition at 10 K was suppressed evidence for a gradual freezing 
of the orbital degrees of freedom has been observed 19 . However, in 
these crystals the conductivity was significantly enhanced compared 
to the samples investigated here and Maxwell- Wagner like relaxa- 
tions extended down to the lowest temperatures 41 . Measurements of 
the dielectric constant as a function of temperature and magnetic 
field had not been reported. Recently the occurrence of multifer- 
roicity in the orbitally ordered phase of FeCr 2 S 4 has been dis- 
cussed 42 . Figure 4(a) shows the temperature dependence of the 
dielectric constant s' of FeCr 2 S 4 around the orbital ordering transi- 
tion at Too ~ 10 K for frequencies between 300 Hz and 10 kHz and 
external magnetic fields up to 1 T. In these measurements the 
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magnetic field was oriented parallel to the electric field, H 1 1 E. For all 
frequencies and magnetic fields the dielectric constant decreases on 
decreasing temperatures reaching constant low-temperature values 
below 6 K in the orbitally ordered state. 

The orbital ordering transition in the 10 kHz-curve is indicated 
by a kink in the dielectric constant and is present for all applied 
magnetic fields. At lower frequencies the transition is masked by 
the increase of the dielectric constant towards higher tempera- 
tures due to an extrinsic Maxwell-Wagner relaxation 41 . The 
strong increase of the dielectric constant in external magnetic 
fields of 1 T at low temperatures indicates a magnetic-field 
induced effect of the polarization and will discussed further in 
the following. 

In Fig. 4(b) we plot the change of the dielectric constant 
Ae(H) =e'(H) — e'(0) induced by an applied magnetic field for the 
configurations H j | E and H _L E. The behavior is clearly nonlinear 
for both configurations. Notably, a positive contribution occurs for H 
1 1 E and a smaller negative one for H _L E. In both configurations we 
observe weak hysteretical behavior in low magnetic fields (not 
shown) and a change of slope in the applied magnetic field range 
4.5 T < H < 5.5 T (hatched area in Fig. 4(b)) were observed. The 
latter effect is also documented in the neutron scattering and mag- 
netostriction results discussed above. For H || E we plot As(ff) 
together with the magnetization M(H) of the sample in Fig. 4(c). 
Both curves can be scaled to fall on top of each other indicating 
Ae(ff)ocM, which even holds in the field range 4.5 T < H < 
5.5 T. This relation is confirmed in Fig. 4(c), where As(H) is plotted 
directly vs. M(H) and reveals a linear correlation between the two 
quantities with a slope of 4.1 (f.u.//.( B ). 

To understand this magnetic-field dependence we assume that the 
polarization of the sample as a function of electric and magnetic field 
may be written as 



where P°(E, H) is a possible spontaneous polarization, which may 
depend on the applied fields, denotes the electric susceptibility 
tensor and a the tensor of the linear magnetoelectric effect. Higher 
order terms in E and H are neglected. As the polarization exhibits a 
strongly nonlinear dependence on the applied field, we also neglect 
the linear magnetoelectric effect in the following. 

In the dielectric measurement the dielectric tensor & r is evaluated 
by assuming 

D = E 0 E + P = E 0 £rE. (2) 

with P = P°(E, H) + Eo^TL The presence of a static electric polar- 
ization P°(E, H) may then result in an effectively magnetic-field 
dependent £ r . Here, investigating a polycrystalline sample and E = 
(0,0, E z ) = E z e z along the z-direction in the laboratory system we can 
distinguish two cases for the dielectric constant according to the 
relative orientation of the applied magnetic and electric fields: 

Ae,| (H) = 8 0 " 1 (P° (E Z ,H Z ) - P° (£ z ,0) -e z ) /E z (3) 



Ae_l (H) = 8" 1 (P° (E z ,H y ) - P° (E z ,0) -e z ) /E z 



(4) 



P(E,H) = P° (E,H) + e 0 /E + aH - 



(1) 



The experiment shows that Aej| (if) > 0 while Ae± (H) < 0, which ori- 
ginates from a change of the polarization component along the z-axis 
in the presence of the applied magnetic field. Applying the magnetic 
field perpendicular to the electric field obviously reduces the projection 
of the total polarization on the z-axis and leads to Aei (H) < 0. 

Given the linear relation between Aej (H) and the magnetization 
M(H) we conclude that within our assumptions the z component 
P°(£ z , H z ) • e z of the polarization and the corresponding ferroelectric 
domains undergo the same changes as the magnetic domains and 
that ferroelectric and ferrimagnetic domains are linked to each other 
in FeCr 2 S 4 . 

A similar situation has recently been reported in single crystals of 
the multiferroic spinel CoCr 2 0 4 32 ' 43 , which is comparable to FeCr 2 S 4 : 
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Figure 4 | (a) Temperature dependent dielectric constant e' of polycrystalline FeCr 2 S 4 for selected frequencies. The measurement were performed in 
zero-magnetic field and with applied external magnetic fields H (0.05 T and 1 T) parallel to the direction of the electric field E. (b) Magnetic field 
dependence of the polarization contribution AP(H) = P(H) — P(0) induced by the applied magnetic field (measured at v = 1 kHz and T= 5 K) for H || E 
and H ± E. (c) Direct comparison of AP vs. H for H || £and magnetization Mvs. H. (d) Plot of AP(ET) vs. M(H) showing a linear correlation of the two 
quantities. 
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it is isostructural at room temperature, becomes a collinear ferri- 
magnet below about 95 K and ferroelectric below 27 K as a result 
of a conical-spiral magnetic ordering 32,43 . At the latter phase trans- 
ition the dielectric constant in CoCr 2 0 4 reportedly exhibits a small 
cusp similar to the effect shown in Fig. 4(a). For CoCr 2 0 4 a spin- 
driven origin has been proposed as the origin of multiferroicity. 

Although the observed field dependence of the dielectric constant 
is not a direct proof of the existence of a static polarization compon- 
ent and multiferroicity, we believe that in FeCr 2 S 4 a similar mech- 
anism might be active. In our opinion, the clearly linear correlation 
between Ae(H) ccM(H) shown in Fig. 4(c) strongly points towards a 
multiferroic ground state of FeCr 2 S 4 . However, the actual spin 
arrangement of FeCr 2 S 4 below 50 K remains unsolved. Assuming a 
scenario where ferroelectric and ferrimagnetic domains are clamped, 
the field dependence of ASj (H) in Fig. 4(c) in fields below 1 T can be 
understood in terms of ferroelectric domain reorientation and align- 
ment effects, which are induced by the alignment of magnetic 
domains. The change of slope at 4.5 T < if < 5.5 T is then related 
to the changes in the magnetization and the increase of the magnetic 
moments of the Fe- and Cr-ions assigned to the overcoming of the 
magnetocrystalline anisotropy. 

Summary and conclusions 

Using neutron diffraction we observed a reduction of the cubic lattice 
parameter below T 00 ~ 10 K in agreement with bulk measurements 
in thermal expansion. This transition was previously attributed to an 
onset of orbital order 12 and is not altered under an external magnetic 
field as demonstrated by our neutron diffraction and specific heat 
measurements. In the field range 4.5 T < H < 5.5 T an increase of 
the magnetic moments of the Fe- and Cr-ions appears. The bulk 
magnetization accordingly exhibits a change of slope, and the mag- 
netostriction causes an expansion of the sample along the direction 
of the applied field. These changes are attributed to the magnetocrys- 
talline anisotropy, which is overruled by the external magnetic field, 
inducing an isotropic ferrimagnetic state for fields larger than 5.5 T. 
In this isotropic regime, the magnetic moments of the Fe and Cr 
reach the expected spin-only values, suggesting that transverse com- 



15 



12 



H 
X 

o 

a. 



FeCr 2 S 4 



00/ 

isotropic FiM/ 
ferroelectric 



-» — O 



OOI non-collinear FiM/ 
ferroelectric 



ou 

non-collinear 
FiM 




10 



100 



Temperature (K) 



Figure 5 | H — T phase diagram of FeCr 2 S 4 (logarithmic temperature 
scale). OO: orbital order (as determined in Ref. 12), FiM: ferrimagnet, OL: 
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magnetostriction AL(H, T) (spheres) and the specific heat C P ( T, H) (stars) 
and are in accordance with our neutron diffraction data presented in Fig. 2. 
Values from the magnetization M( T, H) (diamonds) were taken from Ref. 
31. Lines are to guide the eyes. 



ponents of the magnetization are suppressed together with the mag- 
netocrystalline anisotropy. 

The dielectric measurements revealed a strong dependence on the 
applied magnetic fields. In case of parallel applied electric and mag- 
netic fields, a linear correlation between the dielectric constant and 
the magnetization measured in magnetic field is established. We 
interpret this behavior as an indication for the existence of a ferro- 
electric polarization in FeCr 2 S 4 , where ferroelectric and magnetic 
domains are strongly coupled to each other. To corroborate such 
an interpretation we suggest that future studies may be directed to 
search for structural features of a polar symmetry group in FeCr 2 S 4 . 

Based on these results, we show a modified H — T phase diagram 
of FeCr 2 S 4 in Fig. 5. In the orbitally ordered state below about 10 K, 
we assign the observed features to a ferroelectric ground state and, 
hence, regarded FeCr 2 S 4 as a multiferroic material like the related 
spinel CoCr 2 0 4 32,43 . Most multiferroic materials are antiferromag- 
netic and the spin order cannot be addressed directly by an external 
magnetic field. This problem would be largely overcome by the pre- 
dominantly ferrimagnetic structure of FeCr 2 S 4 and would hence 
open new routes for basic research and applications in this poten- 
tially new multiferroic material. 

The changes in the field range 4.5 T < H < 5.5 T and for the 
lowest temperatures are assigned to a transition into an isotropic 
ferrimagnetic state. The corresponding necessary magnetic fields 
increase when approaching the orbital ordering temperature from 
below and reach up to about 10 T. Due to the increasing conductivity 
with temperature, no such clear anomalies were observed in the 
orbital liquid regime above 10 K, but it can not be excluded that 
ferroelectric domains, which are only weakly clamped to the mag- 
netic domains, are already present in the orbital liquid regime. 

Methods 

Polycrystalline samples of FeCr 2 S 4 for neutron experiments were prepared by solid-state 
reaction using high-purity elements. Post growth annealing under a chalcogen atmo- 
sphere was performed to achieve a pure stoichiometric ratio, which has been shown 
previously to exhibit the strongest low temperature anomaly 12 . For preparation of dense 
samples for thermal expansion and specific heat measurement, we used a spark-plasma- 
sintering (SPS) technique as described in Ref 12. While both samples exhibit the orbital 
ordering transition at about 10 K, the sample obtained by SPS showed a reduced 
ferrimagnetic transition temperature of 165 K in comparison to the one used for the 
neutron diffraction study with T c — 180 K. Thermal-expansion, magnetostriction and 
dielectric experiments were carried out in a home-built setup using a capacitance bridge 
(AH 2700 Hagerling) for temperatures from 1.5-300 K and in external magnetic fields 
up to 14 T. The sample had the shape of a platelet and the expansion was measured 
along the direction of the applied magnetic field, which is perpendicular to the platelet. 
Dielectric measurements were perfomed using silver-paint contacts applied on both 
sides of the platelike sample, the magnetic field was applied perpendicular and parallel 
to the electric field. The heat capacity was measured in a physical properties mea- 
surements system (PPMS, Quantum Design) for temperatures between 1.8 K and room 
temperature and in external magnetic fields up to 9 T. 

High resolution powder neutron diffraction measurements were performed using 
the instrument Echidna, located at the Bragg Institute, ANSTO, Sydney, Australia. 
Echidna is equipped with 128 3 He linear position sensitive detectors that are scanned 
in position to produce high resolution diffraction patterns with an angular step of 
0.05°. A Ge monochromator was aligned on the (331) reflection with a take-off angle 
of 140 to select a wavelength of 2.4395 A with a calculated accessible Q-range of 0.2- 
5.1 A -1 . Collimation ensured a minimum FWHM resolution of ~0.4°. To achieve a 
good signal to noise ratio, each diffraction pattern was collected for 6 hours. The 
symmetry of the FeCr 2 S 4 crystal and magnetic structures were characterized by 
performing a least-squares Rietveld refinement of the powder neutron diffraction 
data, using the fullprof software suite 44 . 

The high-intensity powder neutron diffractometer Wombat, located at the Bragg 
Institute, was used to conduct temperature and applied magnetic field dependency 
studies of the magnetic and crystal structures of FeCr 2 S 4 . Magnetic fields up to 7.5 T 
were applied using a vertical closed-cycle cryomagnet (Cryogenic Ltd.). Wombat uses 
a monolithic position-sensitive 3 He detector, enabling fast collection times with little 
penalty to resolution. A vertically focused Ge monochromator was set to align on the 
(113) reflection, with a 90° take-off angle, giving an incident wavelength of 2.4 1 A and 
an accessible Q-range of 0.5-4.9 A -1 . No collimation was used to maximize the count 
rate, resulting in a required 10 minute collection time per scan. 
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